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Abstract 

The (f>(k±) ~ kj_ relation is called polarization structure. By density functional calculations, we 
study the polarization structure in ferroelectric perovskite PbTi03, revealing (1) the k± point that 
contributes most to the electronic polarization, (2) the magnitude of bandwidth, and (3) subtle 
curvature of polarization dispersion. We also investigate how polarization structure in PbTi03 
is modified by compressive inplane strains. The bandwidth of polarization dispersion in PbTi03 
is shown to exhibit an unusual decline, though the total polarization is enhanced. As another 
outcome of this study, we formulate an analytical scheme for the purpose of identifying what 
determine the polarization structure at arbitrary k± points by means of Wannier functions. We 
find that <f>(k±) is determined by two competing factors: one is the overlaps between neighboring 
Wannier functions within the plane perpendicular to the polarization direction, and the other is the 
localization length parallel to the polarization direction. Inplane strain increases the former while 
decreases the latter, causing interesting non-monotonous effects on polarization structure. Finally, 
polarization dispersion in another paradigm ferroelectric BaTiC>3 is discussed and compared with 
that of PbTi0 3 . 

PACS numbers: 77.22.Ej, 77.80.-e 
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I. INTRODUCTION 



Electric polarization is a key quantity for computing and understanding technologically- 
relevant effective charges, dielectric and piezoelectric responses that are the derivatives of 
polarization with respect to atomic displacement, electric field, and strain, respectively. []J 
Polarization also plays an important role in the methodology development of the theory 
dealing with finite electric fields in infinite solids, by minimization of the free energy F = 
U - E ■ P Total electric polarization consists of electronic contribution (P e i) 

and ionic component (Pi on ) ■ Computing the latter component is straightforward using point 
charges, while calculating the electronic polarization is not. Today P e i is calculated using the 
sophisticated modern theory of polarization jf], Qj. According to the theory, P el corresponds 
to a geometrical phase of the valence electron states, 

2 f 

Pel = 7^1 J dk ± <f)(k ± ) , (1) 

where 

M r G \\ d 

0(*±)=<X)y o dk wKk\^Ku) (2) 

is the Berry phase of occupied Bloch wave functions u n ^. Subscripts || and _L mean parallel 
and perpendicular to the polarization direction, respectively. Practically, to carry out the 
P el calculations, the integral in Eq.([T]) is replaced by a weighted summation of the phases 
at sampled discrete Appoints (Monkhorst-Pack scheme 8[], for example) in the 2D k±_ plane, 
namely, P e \ = Ylk x UJ (k±)(f)(k±) with weight Ylk ± UJ {k±) = 1- The polarization at individual 
k±, 4>(k±), is calculated as the phase of the determinant formed by valence states at two 



neighboring ks on the feii string as 



j-i 

(j>(k ± ) = Im{ln JJ det((u fci>m |u fcj+lin ))} . (3) 

Defined as such, the total polarization P = P ion + P ei could be uniquely determined and 
gauge independent up to a modula constant. In Eq.([T]) one sees that, it is the (f>(k±) phases 
at different k± that determine the electronic polarization. The purpose of this work is to 
study the properties of </>(k±). 

The physical significance of the <f)(k±) quantity can be understood by analogy. It is well 
known that band structure, which describes the relation between single-particle orbital en- 
ergy and electron wave vector k, is very useful for understanding electronic, photo-excitation, 
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and photoemission properties in solids[9|. The 4>(k±) ~ k± relation may be similarly termed 
as "polarization structure" , or "polarization-phase structure" . Electron states in band struc- 
ture can be changed by photo-excitation or emission. The fcj_-point polarization phase can 
be altered by electric fields, which act as a possible excitation source for electrical polar- 
ization. Note that electrical fields do not alter the electron wave vector (k±) perpendicular 
to the direction of the field, and thus k± remains a conserved quantity. The field-induced 
variation of (j)(k±) in fact manifests the fc^-dependent polarization current. As a result, the 
relevance of polarization structure to electronic polarization is like the band structure to 
electronic properties. 

Furthermore, understanding the 4>{k±) quantity is of useful value from both fundamental 
and computational points of view. Fundamentally, this quantity is determined by the Bloch 
wave functions, not in the ordinary sense of spatial distribution, but through the interesting 
aspects of the Berry's phase of occupied manifold of electron states. Studying how </>(fc_i_) 
depends on k± may yield better understanding of electron states, as well as the rather 
intriguing connection between these states and their contributions to polarization in insulator 
solids. Computationally, we first recognize that the (j>(k±) phase computed from Eq.([3]) 
always produces a value within the principal range [0, 27r]. In reality, depending on the 
dispersion of 4>(k±) as a function of k±, it is possible that the phases for different k points 
fall in different branches. In other words, the true <p{k±) values may fall in the principal 
range for some k± points (Let us denote this set of k± points as kf ), while falling out of 
the principal range for other k± (to be denoted as k^). We find numerically that this 
indeed happens for real materials particularly when polarization is large; a specific example 
is given in section II. When this occurs, one must not artificially shift the phases of the 
k^ into the principal range, as computers do according to Eq.(j3J). Though this shift makes 
no difference to the polarization phase of individual k± points, it will alter the total P e i 
polarization, yielding spurious magnitude of polarization. Only when the phase of every k± 
is shifted by a constant 2tt will the total P e \ polarization remain equivalent. To find out 
which kj_ may generate a phase not in the principal range, one in principle should compute 
the whole dispersion structure of polarization and then map out the 4>{k±) for all k± points 
based on the assumption that the (f>(k±) phase is a continuous function of wave vector k±, 
which makes it important to study the properties of the <p(k±) phase as a function of k±. 

Despite the relevance, the dispersion structure of polarization is nevertheless not com- 
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pletely understood. More specifically, (1) little is known about what determine the 4>(k±) 
phase at individual k±. In Eq.([3]), <p{k±) is determined by the wave functions of a string of 
fcii points, not just a single k. As a result, the answer to the question is highly non-trivial. 
(2) For a given ferroelectric substance (say, the prototypical PbTiOs), it is not clear which 
k± exhibits the largest polarization contribution. Does the T point always contribute most 
or least? (3) We do not know if the Berry's phases at different kj_s share a similar value 
or are very different from each other, that is, a problem concerning the dispersion width of 
the polarization structure. Slightly more intriguing, one may wonder along which direction 
the ~ k± curve shows the largest dispersion? (4) Even for two commonly studied ferro- 
electrics, BaTiOs and PbTiOs, we do not know how different or similar their polarization 
structures are. 



Recently, there is another active field in the study of polarization, w 



rich concerns the 
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131 ] . This tunability 



use of inplane strain to tune the ferroelectric polarization 1(3, 
stems from the fundamental interest in the strain-polarization coupling. Imposed under 
inplane strain ferroelectrics subject to modifications of chemical bonds and/or charge trans- 
fer, thereby the interaction between atoms is altered. It has been known that a compressive 
inplane strain tends to enhance the total polarization. But the amplitude of enhancement 



was found to be highly material dependent. 
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13] Considering the importance of the strain 



effects, one might want to know how the 4>(k±) phase from each k± can be influenced by 
strain. Strain effects on the polarization dispersion remain largely unknown, however. It 
would be of interest to examine how the strain may tune and modify the dispersion of polar- 
ization structure. Specific questions on this aspect are: in what manner would the inplane 
strain change the relative contributions and curvatures at different k±, and how the band 
width of the dispersion curve is to be altered. 

With these questions in mind, we here study the dispersion structure of the polarization 
in ferroelectric perovskites, as well as its dependence on inplane strains. Two complemen- 
tary approaches (first-principles density functional calculations and analytical formulations) 
are used. By means of analytical formulation, we aim at a better understanding on what 
specific quantities and/or interactions determine the polarization at individual kj_ point. 
Our calculations reveal some useful knowledge on the polarization structure in perovskite 
ferroelectrics. For example, the largest 4>{k±) contribution is shown not to come from the 
zone center, but from the zone boundary. We also find that the polarization curve in PbTiOs 
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is notably flat along the T — X\ direction, and exhibits, however, a strong dispersion along 
the T — X2 axis. Our theoretical analysis further reveals that the flat dispersion along the 
r — Xi direction is caused by a small amount of participation from the nearest-neighbor in- 
teraction between the Wannier functions. Finally, the present study also demonstrates some 
rather interesting differences in PbTi03 and BaTiOa, in terms of the polarization structures 
as well as their strain dependences. 



II. THE POLARIZATION STRUCTURE OF LEAD TITANATE 



We first present the density functional calculations on the polarization structure in 
PbTiOa. In its ferroelectric phase PbTi03 is tetragonal (|ai| = |a.2 1 = a, \a&\ = c) and 
possesses a large spontaneous polarization. The polarization is along the c-axis direction, 
perpendicular to the k± plane. Calculations are performed within the local density approxi- 
mation (LDA) [14]. We use pseudopotential method with mixed basis set [15] . The Troullier- 
Martins type of pseudopotentials are employed [16| . Details for generating pseudopotentials, 
including atomic configurations, pseudo/all-electron matching radii, and accuracy checking, 



were described elsewhere 



ITT ] . The energy cutoff is 100 Ryd, which is sufficient for conver- 



gence. The calculations are performed in two steps: the optimized cell structure and atomic 
positions are first determined by minimizing the total energy and Hellmann-Feynman forces, 
and after the structural optimization, the polarization dispersion of (j)(k±) is calculated us- 
ing the modern theory of polarization. (j, Q] Our LDA-calculated inplane lattice constant 
for unstrained PT is a=3.88A, with c/a — 1.04, both agreeing well with other existing 
calculations. 

Figure Ufa) shows the reduced 2D Brillouin zone that the k± points sample over. The 
calculated (f> phases at individual k± points along the r — > X\ — * X2 — > T path are given 
in Fig{T](b). Reciprocal-space coordinates of X\ and X2 are k± = (jr/a,0) and (it /a, it /a), 
respectively. The dispersion curve is rigidly shifted such that the phase at T is taken as the 
zero reference. 

Before we discuss the specific results in FigJU we need to point out that the shape of this 
&j_-dependent phase curve is translation invariant. As is known, the electronic polarization 
alone can be an arbitrary value, if the solid is uniformly translated with respect to a fixed 
origin of coordinates. Though different translations will change the absolute location of the 
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polarization-dispersion curve, the shape of the curve remains unaffected, however. This can 
be easily illustrated by analyzing the change in the <f)(k±) phase when one displaces the 
solid arbitrarily. Let the wave function of the original system be if> z(r) = e u Ar), where 
u nk(f) = u nk( r * Jr ^) • Now, we displace the solid by an arbitrary vector while the origin 
of coordinates is fixed. Let us denote the original system using script A and the displaced 
system using script B, so rg = ?a + Tq. The wave functions of the displaced system satisfy 

V£fo) = = 1&fB - ro) • (4) 

Thus we have u B Jr B ) = e~ lkr °u A JrB — ro)- Substituting this relation into Eq.([2J) or Eq.Q, 
one can obtain that the 4>(k±) of the displaced system is 

<P B (k ± ) = <j ) A (k ± )+f .G ll N^ d , (5) 

where N^ d = M is the number of bands occupied by electrons. The phase differences 
between the A and B systems are thus a constant, independent of k±. 

Several observations are ready in Fig{T](b): (1) The largest (p(k±) polarization does not 
come from the zone-center T point. Rather surprisingly, the largest 4>{k±) phase is from the 
X% point which lies at the far end of the BZ. (2) The polarization curve is flat along the 
r — X\ line, showing only a small dispersion. On the other hand, the dispersion becomes 
very large along the V — X 2 direction. (3) At k± points of high symmetry (such as T, A 1; 
or X2), the curve in FigfjTb) has zero slope, similar to the electron band structure. (4) The 
dispersion of polarization also shows subtle details which could not be easily understood. 
For example, there is a local (though not very pronounced) maximum along the Ai line, 
making the X\ point a local minimum in both r — X\ and Xi — X 2 directions. 

Our calculations further reveal that, despite the fact that the polarization in Figfjjb) 
exhibits substantial kj_ dependency, the dispersion width (~0.6) is much smaller than 2n. 
This finding is important for the following reason. As described in the introduction, if 
the differences of the (f>(k±) phases at different k± points are greater than 2tc, one would 
encounter a difficulty in determining which branch of phase a specific k± point should be 
assigned. This difficulty can be avoided only after the phases of all k± points are mapped 
out. Fortunately, the result in Figfjjb) tells us that the phase contributions from different k± 
points are fairly close, and the differences are far less than the critical value of 2tt that may 
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cause the above difficulty. Nevertheless, we should point out that even a small polarization 
dispersion as in FigfjTb) may still give rise to spurious results on total polarization. To 
illustrate this, we displace all five atoms in PbTiC>3 along the polar c-axis by a distance 
Zq. FigfSJ^a) shows the total (electronic +ionic) polarization, computed from the geometric 
phase, as a function of the displacement z (in unit of c). Intuition tells us that the total 
polarization should be uniquely determined and translationally invariant. However, we see 
in Figj2](a) that unphysical discontinuity happens for some zq points, and this discontinuity 
shows up periodically. To understand what causes the discontinuity, we examine the phase 
contributions from individual k± (sampled according to the Monkhorst-Pack scheme Q]), 
depicted in FigfS^b). Figure 0(b) shows that the individual- k± phases indeed are a periodic 
function of Zq, explaining why the discontinuity in Figf5](a) is periodic. Here it may be useful 
to comment briefly on the length of the periodicity. One might think that by displacing the 
solid by a distance of c in the c-axis direction, the 4>{k±) phase would change by a value of 
2tt. However, the periodicity in Figj2]is much smaller than c. The explanation is simple. As 
a matter of fact, in real space the individual 4>(k±) has a periodicity of N } cc c (instead of c), 

band 

which for PbTi03 the periodicity is 0.0455c because Ng£° d = 22. This is indeed consistent 
with the numerical calculation in PT (FigfSb). The length of periodicity can be seen from 
Eq.Q, showing that, whenever r = N ™ cc R\\ (n is an arbitrary integer and R\\ is the lattice 

band 

vector along the G\\ direction), the (p B (k±) and (fi A (k_i_) differ by (f> B (k±) = (p A (kj_) + 2%n. 
FigEJ^b) also reveals the reason responsible for the discontinuity of the total polarization. 
Spurious discontinuity occurs when the 4>(k±) phases of some (but not all) individual k± 
exceed 2n [Figfj^b)]. Under this situation, computers incorrectly shift the phases of these 
k± points back to the principle range, yielding spurious total polarization. According to our 
experience, spurious polarization often takes place in two circumstances: one is for materials 
of very large polarization, such as tetragonal BiScC>3, and another is when atoms in the unit 
cell are translationally shifted. Given the small bandwidth of the <fi(k±) dispersion, it is 
now straightforward that, by using different f*os, we can avoid the spurious polarization. 
However for some materials, if the dispersion width from different k± points is larger than 
27r, one may have to rely on the continuity of the 4>(k±) phases, and map out the phases of 
individual k± points over the whole two-dimensional k± plane in order to find the correct 
phase branch. 
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III. STRAIN DEPENDENCE OF POLARIZATION STRUCTURE 



An important property of ferroelectrics is that the polarization is strongly dependent 
on strain. While strain can change the total polarization, response of the polarization 
dispersion structure to strain could also be an interesting problem. Here we investigate the 
response of the polarization structure under inplane strain in PbTiC>3. For each in-plane (a) 
lattice constant, the out-of-plane c lattice constant and atomic positions are fully relaxed, 
by minimizing the DFT total energy. The polarization structure is then determined using 
the optimal structure. 

Figure [3] shows the phase dispersion curves for PbTiOa at different inplane lattice con- 
stants. All curves are shifted so that the phase at T point is zero, in order to conduct 
direct comparison. Three conclusions can be drawn from Figj3j (1) The relative phase, 
— changes drastically for X 2 , but not so significantly for X\. (2) At increasing 

strain, (or smaller inplane a constant), the bandwidth of the dispersion initially changes 
very little when a = 3.84A, and then starts to decrease upon further increasing strain to 
a = 3.80A. The decline of the dispersion bandwidth is rather surprising, since a compres- 
sive inplane strain is known to enhance the total polarization in PT. The decline is also 
counterintuitive when one considers that the decreasing inplane lattice constant makes the 
atom-atom coupling stronger within the inplane directions, and should therefore have in- 
creased the bandwidth. One possible reason that may cause the decrease of the bandwidth 
is given in the next section. As a result of the declining dispersion, the polarization curve 
becomes notably "flat" at small a = 3.65A. (3) The curvature of the dispersion also shows 
subtle changes, featured by the fact that a new dispersion minimum appears along the X 2 — T 
line at large strain. As a consequence, the dispersion curvature [i.e., the second derivative 
Vl P(k±)] & t r point alters its sign from being positive (at large a) to negative (at small 
a). Furthermore, the local maximum between r — X\ for unstrained PT turns into a new 
minimum at large inplane strains. Meanwhile, the X\ point changes from a minimum into 
a saddle point, when strain increases. 

The calculations thus reveal that, while inplane strain has been previously known to 



12 1 and sometimes re- 



introduce interesting modifications (sometimes markedly enlarged 
markably small [13j) to the total c-axis polarization, its effects on the polarization dispersion 
at individual k± points appear to be even richer, showing that the polarization structure 
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indeed worths studying. The subtle response of the polarization structure, as predicted 
above, indicate that there is new and rather complex physics behind the results in Figj3l 
While we know that the strain-induced changes in the polarization dispersion must be asso- 
ciated with the fundamental modification of electron wave functions, we also have to admit 
that the DFT results obtained in our numerical calculations are puzzling, and an intuitive 
understanding of the results is difficult for two reasons. First, this is an early attempt to 
investigate the polarization structure, and there is not much previous understanding in the 
literature. Second, although Eq.pj) and Eq.Q allow us to compute precisely the polariza- 
tion of individual k±, a direct and more intuitive connection between <f)(k±) and Bloch wave 
functions is hard to capture from these equations. As a result, it would be very helpful if one 
could find an alternative way to understand the polarization structure and the computation 
results. For instance, what determines the polarization at individual k± point, and why 
<f)(k±) maximizes at the X 2 point? In the next section, we attempt a scheme which we wish 
to be able to offer a more intuitive understanding of the polarization structure. 

IV. WANNIER FUNCTION FORMULATION OF POLARIZATION STRUC- 
TURE 

As mentioned above, Eq.(T2]) and Eq.(J3]) give us little intuitive sense on the direct k± 
dependence of the Berry's phase. In order to get more insight, we use Wannier functions to 
analyze the polarization structure. Previously, Wannier functions have been found very use- 
ful in analyzing real-space local polarization^, Q]- Here we employ the Wannier-function 
approach for a different purpose, namely to understand the fc^-dependence of the polariza- 
tion structure. The Wannier functions are defined as 





or 




(7) 



where R runs over the whole real-space lattice vectors. By substituting Eq.([7]) into Eq.(T5]) 
and carrying out analytically the integral over k», it is straightforward to derive, for tetrag- 
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onal perovskites, the polarization at individual k± as 

?7T - f 

<j>(k ± ) = —J2J2 W:(f)W n {r-R ± y k ^dr (8) 

° R ± n=l J 

where fji is the projection of vector r along the polarization direction, R± is the projection of 
lattice vector R onto the plane perpendicular to the polarization direction. For convenience 
of discussion, we separate the sum over R±_ into the R±_ = term and the rests, 

9 M r 

<j>$i_) = fo + —J212 nW:{?)W n {f - R ± )e lk ^dr (9) 

R ± ^0 n=l 

where for R± = 0, (f> = ^^ = x J (^\\W*(r)W n (r)df is the phase contribution from the 
same unit cell. Eq.([9]) is the basis for understanding the polarization structure. From this 
equation, we observe the following. 

First, it is now clear that the /c^-dependent part of (p(k±) comes only from the R± ^ 
terms, which correspond to the overlap of the Wannier functions in neighboring cells. In 
other words, the k± dependence of the (j>(k±) phase results from the overlap of the Wannier 
functions of different cells that are displaced by R± from each other within the plane that 
is perpendicular to the direction of polarization. While the choice of the Wannier function 
is known to be non-unique due to the gauge uncertainty, the sum of the Wannier-function 
overlap over occupied bands is a uniquely defined quality which does not depend on the 
gauge. It is this quantity that determines the shape of the polarization structure. 

Second, Eq.Q explains why the bandwidth of polarization dispersion is often much 
smaller than 2n. Since only the second term in this equation is k± dependent, and since the 
Wannier functions are generally well localized compared to the size of unit cell, one expects 
the overlap W*(r)W n (r — R±) to be much smaller than unity for R± ^ 0. This is consistent 
with our numerical results in FigFTJ, namely, <f>(k±) — <po ~ 0.6 2n. 

Third, since the dispersion in <p(kj_) comes from the overlap of the Wannier functions 
between cells of different R±s in the x?/-inplane directions, it explains why the polarization 
structure is very sensitive to inplane strain, where by changing inplane lattice constant, the 
distances between neighboring cells are effectively altered. Meanwhile, we recognize that a 
precise understanding of how the bandwidth depends on the inplane strain is not as simple 
as one might think. Naively one tends to think that, with the decline of inplane lattice 
constant, the dispersion is to increase, since the overlap W*(r)W n (f — Rj_) increases when 
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R± decreases. This will lead to the widening of the polarization dispersion width, which is 
opposite to what we found in Figj3j This puzzling contradiction can be resolved by noticing 
that, in addition to being dependent on the overlap strength between W n (r) and W n (f— R±) 
within the perpendicular plane, the dispersion width also hinges on the localization length 
(/jj VF ) of the Wannier functions along the direction parallel to the polarization, as a result 
of the f\\ operator in Eq.([9]). With the increasing inplane strain, the is to shrink. We 
thus see that the bandwidth of polarization is determined by the balance of two competing 
factors between the increasing Wannier-function overlap and the decreasing Zj^ localization 
length. When the latter dominates, the bandwidth declines as we have seen in Figj3] from 
numerical calculations. 



V. CURVE ANALYSIS 

With the general understanding of the polarization structure in the above section, we 
next attempt to determine analytically the polarization dispersion specifically for PbTiC>3, 
aimed to obtain further insight into the important details of the polarization structure. As 
will become clear later, our analysis in the following also explains what determines the <j){k±) 
polarization at special points of T, X\ and X 2 . We begin by defining parameters 

M „ 

r\\W*{r)W n (r- R ± )df , (10) 



t(Ri 



2tt 

~ i 

n=l 



and then, 



(11) 



R± 

For dielectrics of insulating nature, Wannier functions are highly localized, and decay 
exponentially with the distance 2C|, l21[. As a result, t(R±) also decay quickly with the 



increase of \R±\, so we can adopt the tight-binding like approach and consider only several 
R±s that correspond to some nearest neighbors (NN). We consider up to the 2 nd NNs, where 



Ri 



(0 0) on site 

(±a 0) INNs 

(0 ±a) INNs 

k (±a ±a) 2NNs 
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Taking advantage of tetragonal symmetry, we can rewrite Eq. liTTj) as 

(j>(k±) = t + 2ti[cos(fcia) + cos(/c 2 a)l 

(12) 

+2* 2 [cos(/ci + k 2 )a + cos(/ci — k 2 )a] , 

where *j is the i th NNs contribution defined in Eq.(lO), and k± = {kx : k 2 )- This expression 
gives us a more direct sense of the <j>(k±) ~ polarization dispersion, approximated to the 
second nearest neighbors. At special fc^ points of T, X±, and X 2 , the phases are 0(T) = 
t + 4ti + 4* 2 , 4>{Xi) = t - 4* 2 , and 0(X 2 ) = t - 4*i + 4* 2 , respectively We could thus 
clearly see that the *o term, corresponding to R± = 0, acts to rigidly shift the polarization 
curve as a whole. Meanwhile, the phase relative to the T (i.e., the dispersion) is determined 
by the t\ and t 2 quantities, and more specifically, 

<f>{X x ) - 4>{T) = -At, -8t 2 , 

<P(X 2 ) - <f>(T) = -8t a . (13) 

These equations are useful, since they tell us that (1) the relative height at X 2 (which 
contributes most to the polarization in PT), 4>(X 2 ) — (fi(T), is determined by *i, associated 
with the overlap of the Wannier function in the I st NNs. *i <0 for PbTi03 in equilibrium. 
(2) Under the assumption that t 2 is negligible, 4>{X 2 ) —<p(T) will be larger than <p(Xi) — <p(T) 
by a factor of 2. 

Within the second nearest-neighbor approximation, one can further determine analyti- 
cally the dispersion along the r — > X\ —>■ X 2 — > F line in the 2D Brillouin zone as 



0(fcj 



t + 2*i + (2*i + 4* 2 ) cos(fcia), for r -> X x with k 2 = 

*o — 2*i + (2*i — 4* 2 ) cos(A; 2 a), for X\ — > X 2 with ki — it /a 

to + 2* 2 + 4*i cos(fcia) + 2* 2 cos(2fcia), for X 2 — > V with k± — k 2 ■ 



The polarization structure could thus be expressed as a simple combination of cosine func- 
tions. 

To examine whether the second-NN approximation is sufficient, we fit the analytical 
results to the numerical DFT calculations to determine the U (i — 0, 1, 2) parameters. Note 
that only <p(kj_)s at three points (i.e., T, X\ and X 2 ) are fitted. The obtained *j values 
are given in Table [0 These values are then used to determine the whole dispersion curve, 
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TABLE I: The fitting ti and t 2 parameters for PbTi03 at different lattice constants. t§ is not 
shown here since it does not affect dispersion. 

a(A) t x t 2 

3.88 -0.072 0.031 

3.84 -0.072 0.032 

3.80 -0.064 0.031 

3.72 -0.031 0.023 

3.65 -0.010 0.016 

shown in Figfljb) for PbTiC>3 in equilibrium structure of a = 3.88A. We could see that the 
analytical curve agrees well with the DFT result, implying that the 2nd NN approximation 
works. On the other hand, some fine structure of the curve (such as the small local maximum 
along the r — Ad) can not be reproduced, where for a better fitting, approximation beyond 
the 2nd NNs would be necessary. 

From Table [I] one can also see how the tj quantities are influenced by inplane strain. 
ti declines substantially as a decreases below 3. 80 A, while t 2 shows a less dependence on 
inplane strain. This makes sense since, by varying the inplane strain, the main effect lies 
in altering the nearest-neighbor interaction among Wannier functions. For a > 3.80A, \ti\ 
approximately equals 2\t2\, confirming the importance of the nearest neighbor interaction. 
For large strains of a < 3.72A, \ti\ and l^l become comparable, for which it is likely that 
higher orders of NNs are also needed. 

VI. COMPARISON WITH BARIUM TITANATE 

It is of interest to compare the polarization dispersions between BaTiC>3 (BT) and PbTiC>3 
(PT), since these two substances have rather different tetragonality, magnitude of polariza- 
tion, and sizes of A-site atoms. For this purpose, we have studied the polarization structure 
in BT, for which a tetragonal symmetry is enforced so that a direct comparison with PT 
can be made. Following the same procedure as for PT, we optimize the cell structure and 
atomic positions of BT at different inplane lattice constants, and calculate the corresponding 
polarization structures. 
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FigJH displays the polarization structure for BaTiOs at different inplane lattice constants. 
Let us first focus on the dispersion of the equilibrium BaTiOs. The LDA-calculated equi- 
librium inplane lattice constant of BT is a = 3.95A. Apart from similarities to PT (e.g., 
(J) maximizes at X 2 ), our calculations reveal some interesting differences between PT and 
BT under zero strain: (1) The BT dispersion curve has a significantly smaller bandwidth 
(~0.42) than that of PT (~0.57). Since the bandwidth is determined by the difference 
<fi(X 2 ) — 0(r), i.e., by t\, a smaller bandwidth indicates less overlapping Wannier's functions 
between nearest neighbors in BaTiOs, which could be explained by the larger inplane lattice 
constant a for BT at equilibrium. (2) Unlike PT, the polarization in BT is not small at 
X\. This again can be attributed to the large inplane lattice constant in BT, which leads to 
a negligible contribution from the 2nd NNs, i.e., t 2 is small in BT. Indeed, we numerically 
found that t 2 is -0.007 in BT, compared to 0.031 in PT. By Eq. ffTB"]) , 4>(Xi) is about half 
of the <fi{X 2 ) value if t 2 is small, which is indeed born out in FigJH (3)As a consequence 
of observation (2), the dispersions of BT and PT along the r — > X\ are not quite similar. 
There is a local maximum between V — X\ for PT, whereas for BT, no local maximum exists 
and Xi becomes a saddle point. 

Upon strain, BaTiOs and PbTiC>3 exhibit sharp difference in their strain dependence 
of dispersion bandwidth. As we saw previously in Figj3j inplane strain causes the band- 
width declining for PbTiCV However, for BaTiOs, a dramatic enlargement in bandwidth 
occurs, when a decreases from 3.95A to 3.85A. The bandwidth maintains a large value at 
a=3.75A, after which it starts to drop. In BaTiOs the polarization dispersion bandwidth 
thus shows an interesting non-monotonous dependence on inplane strain. This characteris- 
tic non-monotonous dependence strongly supports our conjecture that the two competing 
factors determine the bandwidth, as described above in Section IV. When strain is small in 
BT, the overlapping of Wannier functions located at the nearest neighboring R±s plays a 
dominant role, and the increasing overlap leads to a larger \ti\ and thus larger bandwidth. 
As inplane strain becomes large (a < 3.85A), the atom-atom interaction along the c-axis is 
considerably weakened due to elongated c-lattice length. As a consequence, the shrinking 
Ij^ localization length of Wannier functions along the fy direction takes over and becomes 
dominant, giving rise to the declining bandwidth. This, once again, reveals that the po- 
larization dispersion contains rich information. To make more quantitative comparison, we 
replot in FigJS] the strain dependence of the <j>(k±) phases at X x and X 2 , relative to the V 
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point. Figj5]is of some useful value since it allows us to contrast the fc^-specific polarizations 
in two materials at the same fixed inplane lattice constant. The difference between BT and 
PT is thus not related in a significant sense to atom-atom distance, but largely due to the 
overlap of respective Wannier functions. In FigfSl both <p(Xi) and (^{X-z) are seen to be 
far greater in BaTiC>3 than in PbTiC>3, for a fixed a constant. The greater values of 4>{k±) 
in BT could possibly originate from the fact that the Wannier functions in this material is 
more spreading due to the larger size of Ba atom. 

From the comparison between PT and BT, we could see that the polarization structure 
has some common features for materials with similar structure, and meanwhile, some dis- 
tinctions revealing the identities of materials. The common features allow us to understand 
the polarization structure in general, just as for band structure, most III-V semiconduc- 
tors have direct band gaps. Differences in polarization structure manifest the electron wave 
functions and interatomic interactions on microscopic scale. 

VII. SUMMARY 

Two different approaches are employed to study the polarization structure in perovskite 
ferroelectrics. Numerically we use the density functional total-energy calculations and the 
modern theory of polarization. Analytically we formulate a scheme to describe the k± depen- 
dence of the polarization phase using Wannier functions. By parameterizing the Wannier- 
function overlapping, we further identify the quantities that determine the <f>(k±) phases at 
special kj_ points of interest. Our specific findings are summarized in the following. 

For PbTiC>3 at equilibrium, (i) the <fi{k±) phase maximizes at the Brillouin zone boundary 
of the 2D k± plane, not the zone center, (ii) The polarization structure shows little dispersion 
along the T — Xi line. However, the dispersion is large along the V — X2. (iii) The bandwidth 
of the dispersion curve is far below 2ir. The small dispersion considerably eases the difficulty 
in assigning the correct branch of individual k± phase, but caution still needs to be taken 
when the 4>{k±) phase is approaching 2n. 

Analytically, (iv) the expression, Eq.(j2j), is given as the basis for understanding the 
polarization structure. It also explains why the polarization bandwidth is small compared 
to 2-7T. (v) The polarization phase at individual k± is revealed to depend on the competition of 
two factors, namely the overlapping strength of Wannier functions within the perpendicular 
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R± plane and the localization length lY F of these Wannier functions, (vi) Within the 
2NN approximation, the <f)(Xi) and 4>(X2) values in ferroelectric perovskite are found to be 
0(Zi) - <f>(T) = -Ati - 8t 2 , <f>(X 2 ) - 0(r) = -8*i. If t 2 is negligible, the latter is 2 times 
of the former, (vii) When PbTi03 is under compressive inplane strain, the polarization 
bandwidth is found to decrease, whereas the total polarization increases. The declining 
bandwidth implies that the localization length £y of Wannier functions plays a dominating 
role in PbTiCV 

By comparing BaTiC>3 with PbTiOs, we show (viii) the equilibrium BT exhibits a smaller 
bandwidth of 0.42, as compared to the bandwidth of 0.57 in PT. (ix) (f)(Xx) in BaTi03 is not 
small, unlike PT. The difference comes from the fact that t 2 is negligible in BT, leading to the 
result that <p(Xi) is about half of the value of <fr(X 2 ). But in PT, t 2 can not be neglected, and 
acts to offset the t\ contribution, giving rise to smaller 4>(Xi) and flat dispersion along the 
r — X\ line, (x) As BaTiC>3 is under increasing inplane strains, its polarization bandwidth 
displays a characteristic non-monotonous variation by first increasing dramatically and then 
declining. The finding lends a support to the qualitative understanding that two competing 
factors determine the 4>{k±) phase, (xi) When BaTiC>3 and PbTiC>3 are constrained to the 
same inplane lattice constant, the (f)(Xi) and <f){X 2 ) are shown to be significantly larger in 
BT than in PT, unlike the case when two materials are in equilibrium. 

We conclude by pointing out that there are still many aspects of polarization structure 
we do not yet understand. For example, we have not pursued beyond the 2nd nearest 
neighbors to explain the local maximum between L and X\ in unstrained PT. We also do 
not know the physical significance when 4>(Xi) changes from a local minimum to a saddle 
point as displayed in Figj3] for PbTiC>3 under strains. We believe that further analysis of 
the polarization structure could yield better knowledge on the physics of dielectrics. Like 
band structure of solids, we hope that the polarization structure can provide us a new tool 
of studying ferroelectric materials and properties. 

This work was supported by the Office of Naval Research. 
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FIG. 1: (a) The 2D Brillouin zone for the k± plane; (b) Berry's phase at different k± points for 
PbTiOs at equilibrium (symbols: direct calculation results; curve: analytical results). The 4>(k±) 
phase is in units of radian. 
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FIG. 2: (Color online) (a) Total polarization in strained PbTi03 of inplane lattice constant a = 
3.72A as a function of the uniform displacement zq of five atoms; (b) the ^(fej_) phases at six 
Monhorst-Pack sampling k± points as a function of zq. For each c/N^ d change in Zq, the (f>(k±) 
phases change by 2ir. In (b), the <f>(k±) phase curves are enlarged in the right side of the figure for 
zq between 0.044 and 0.048. 
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FIG. 3: (Color online) The (/> phases of different fcj_-points, for PbTi03 under different inplane 
lattice constants. Symbols are direct calculation results; curves are guides for eyes. 
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FIG. 4: (Color online) Polarization dispersions for BaTiOs at different inplane lattice constants. 
Symbols are direct calculation results; lines are guide for eyes. 
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FIG. 5: Dependencies of the (ft(k±) phases at X\ point (left) and at X2 point (right) as a function 
of inplane lattice constant, for PT and BT. 
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